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ABCTRACT 



In a vex'tical plate or t\ibe csxidenacr, the liquid condensate layer 
(filn) flows downward under the influence of gravity forces. Leaninar 
flow generally exists on the t^jpcr part of tl» surface, but txu-bulent 
flow will exist on the lower parts if the condensation rate is suf- 
ficiently high. If tije VS 4 POT velocities are appreciable, the vapor 
exerts a shear stress oti the liquid -vapor interface, v^lch affect 
the heat transfer coefficients. 



The case of flow, \nder the influence of this vapor shear 

sti'esa, has been analysed by Kiwselt (1). By asawalng the "univeiv^sl 
velocity distribution" E^andtl-Eiicuradse, Seban (10) obtained ana- 
lytical predictions for the txarbulent fila, but neglected the effect 
of vapor velocities. This thesis exteiida the previous work by (ievelcp- 
ieg equatloios which Include the vapor shear stress effect on a turbuln •„ 
fils. Two criteria for transition are proposed. Applieatlexi of these 
criteria perialt the present resialts to be ecmbSxved with those of .\*selt 
to give analytical boat transfer coefficients for the case of lapalner 
flow on the \ipper portion of a plate and. turbulent flow on the lower 
portion. 



The results of the analysis are presented lii the forts of dinscasion- 
less plots. Figures 7 to 12, of heat transfer coefficients versus liquid 
Reynolds Jftjoiber at the bottca of the plate and also versus plate length. 
Dimensionless shear stresses and liquid Prandtl numbers ai*e pa^^ssseters 
in these plots. The analytical results ere shown to agree closely 
with Colburn's correlation equation (9) iJi tbe range Np^ ■ 2 to J. 
Eowever, considerable deviation exists la the very hign or very low H_ 
range. The Praadtl EvBbers covered by this analysis extend from 0.01 
to 10. 



An experimental set-up for investigating the hydrodyaaeilcs of a 
falling liquid film is described. Preliminary results Indicate tiait 



the aneblytlcal equatloae axe able to predict fl2a thlcti...3S qu.<^' ac- 
cta*ately f<K* low F.eyuolda flvBibere. However, at falser fteTnold/.. ./ucibera, 
saae deviation was observed. It is co:iclud«i that the Prandtl-iJllniradse 
velocity distribution is generally ade<iu&te for the present purpose 
of predicting heat transfer coefficients, but that ftirther sttrfy is 
necessary to verify this point. 

A procedia*e utilizing a hot wire anetaaaeter is described, ^diich 
appears to hold premise in deteraiaing transition points in the liquid 
film. 
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PAFT I 



83CTI0S A - imoDucnoi^ 



The influ*Dce of a vapor Ghearing stress oa coadenslng fila heat 
transfer coefficients has been ehovn by several investigators (1, 2, 

5, to bo qpilte slgrt if leant. Hovever, this variable is often neglected 
Ir. the interest of a maunageable theoareticaJl analysis vith the result 
that experlEsentally deteitalned values cf heat transfer coefficients 
have often differed froa predicted values by several fold. The dcr/ia- 
tioa has been Euch more noticeable ia the case cf vertical tube conden- 
sers where the condensate flows in film forn on the inside of the tubes. 
This sems reasonfitole when one considers the fact that hig^ vapor vel- 
ocities are «ore easily attainable in this case than in the case of 
horizontal tube condensers or vertical condensers vith the condensate 
in tJiS shell. 

Althou^ higher film coraiensate coefficients can be obtained with 
high vapor velocities it nsust be reambered that for vertical tube con- 
densers, the price of this increase in heat transfer ability ia en added 
pressure drop in the ccsaleneing section. In tlie case of steam conden« 
sera it is also to be noted that the cooling water side coefficient Is 
casually controlling and little can be gained in overall heat transfer 
coefficients by decreasing the condensate resistance to beat flow. Ilow- 
ever, many industrial applications involving the condensation cf organic 
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v«^rs ere coacerned vith condensate coefficients ao high as, or high- 
er than the cooling vxiter coefficients. It is in cases eiich as these 
that reductlcms in size and veight aay be possible throiigh the judicious 
use of this vapor shear effect. 

A condensate film flowing vertically uMer the influence of grav- 
ity and a vapor shear stress lends itself remarkably well to theoret- 
ical attack. It is the p’urpose of this thesis, therefore, to attempt 
to further the previous analyses by deriving expressions for average 
heat transfer coefficients which include both Insinar and turbulent 
cesadensate layers vmder the influence of a vapor shear stress. To sln- 
plify the I'esulting expressions, and at the saiae time give results of 
universal applicability, the case of a flat plate is considered. The 
equations joay also be applied to tubes as Icaig as the thickness of the 
condensate film is small compared to the t\d>e radit^s, and os long as 
the effects of a press'ure gradient may be neglected. 




I 







SECTION B - LW!EJAH FIIM 



In the folloviag analysis tlis ease of a flat vertical plat®. Fig- 
ure 1, is considered. Cmdcnsation is of a ptcre, saturated vapor with 
no noacondenslble gases present. AssiEsiJtions InclMo constant vw.ll 
teaperature and constant flxxld properties. The specific heat, c^, den- 
sity, p, end thennal conductivity, k, may be deteraiined at the arith- 
i^etic average filia temperature. The viscosity, jx, however, should be 
deterained at the special average film temperaature recemended by Col- 
burn (5), and defined by 



' 5A ^ t (1) 

The vapor shear stress, is considered as being directly downward, 

and being constant along the length of the plate. 

The lasiiaar analysis essentially follows that of ffusselt (1), where 
the fllia is in ctreaaline action and is under the influence of a vapor 
shear stress, , as shovaa in Figure 1. 



Wdll 




Vapor 



FIGURE 1 
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B/ ass\sain 5 . steady flow and neglecting lacaaentvja changes, a force bal- 
ance for the shaded element shown gives, 

X - (p ~ Py) - x) - - 0 (2) 

This gives the linmr distribution of shear stress shown in figure 2. 
Since X ae at X « 0, the •wall shear streso is given by 

Xo a (p - p^) x^ + (5) 




X 

FIGTOE 2 

The effect of the Komentua crossing the eletneat faces vlthin the 

liquid Is quite saiall. To include the saoKeatmi effect at the liquid- 

vai)or interface the tena ^xich should be added ’to Equation (2) ia (l/g^)v ^ dC 

o z,x^ 

This represents the z direction velocity of the vapor crossing the li- 
quid -•vapor interface, -which is ■very noerly equal to the liquid surface 
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velocity, v , ead aot th® velocity of the vupor flow in tbo «aia 
vapor streaa. This mosentm tcva. also is usually quite onall. 



Froft ?levton*6 definition of viscosity, wa have for laciinar flow, 

dv 



T 



So 



fhibetitutinE in Equation (2) for L givee. 






a (p . p„)(x„ • X)dX -*■ ^aT(3x 

£ t* *^v'' O ' P 



(^) 



(5) 



This express icsa a»ay he Integrated across the film vith T^^ccmotant, 
giving the velocity profile as 

2i 



^ ^ • ‘"'y> <V ■ ' "5“ *■*' * 



( 6 ) 



The Eass rate of flow per unit width of the vertical pLate, P , 
nay be found by 



r 




p to 



(7) 



Carrying out this integratioa gives P in terms of film thicJcnecs, x , 



r 



-£ 

u 



5 ^2 
, . ®o ^ 

g(p - P,,) — ^ —2^ 



(B) 



Equation (8) may be differentiated with respect to to give a rela- 
tionship between the change in mass flow rate per unit width and the 
change in film thickness, 




e (p - Py) + 




(9) 



A linear taaperatiire gradient may be assumed to exist within the 
condensate film without introducing serious error. Reference (6) gives 
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FIGU^J 3 

If the enthalpy of ^ entering thrcnjgh the upper ®urf«ce of the con- 
trol voltK® is taken at a ssean film temperature. It can be seen that 
the r leaving through the lower surface of the control volme will 
also be at this seme film teaperatiare and hence will not enter into 
an energy balance. It is also evident that the enthalpy change of 
the increment of condcmsate, J ^ , will be equal to the enthalpy dif- 
ference between satxsrated vapor end a liquid at the mean film tempera- 

t 

ture. This may be designated h^^ and has been evalmted by Colburn, 
Millar, and Vestwater (7) os » h^^ + ^/d At. On the basis of 
the previous discussion an energy balance for the control volxsae shown 
in Figure 5 may be inade by equating the heat transfer to the change 
in er.thalpy of the Incremental increase in condensate, iT, 
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(10) 



Since the flte considered is in lasainar motion, the heut trans- 
fer coefficient my be talcen as the thertaal conductivity of the fluid 
divided by the fila thickness, 

» h r>t (il) 

^ ^o 

Eliminating q/A betvecn Eqmtlom (10) and (11) and substituting for 
1 ^ its equivalent from Equation (9), ve have a differential expres- 
sion containing only dz and dx^ as variables, 



k At 



^ p(p - ^ p (12) 



Equation (12) nay be integrated froti the top of the plate where « 0 



when z « 0, to a point x^,z. 



4 



^ ^ p(p . p ) * A. 






(15) 



By rearranging and substituting the Prandtl number, K *» c u/k, the 

Jr^ kr 



resulting expreaslcn is obtained, 

fz ^ (xt-J— = A 



f 1. 5. T- 



N. t (0- " ^ (01, 



The following dimraasianless lengths may be defined as 
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rurther define a dlmensionlees vapor shear stress, , as 









f/-Aj (f) 



‘^'1 



( 17 ) 



Substituting these definitions i:i Equation (l^) gives a isimple acid 

convenient e^cpressloa for the variation of film thickness with length 
♦ 

of the plate, z . 

z* « + V5 (18) 

The usefulness of the dimecsionles© quantities defined by Equa- 
tions (15), (l6), and (17) »ay be more fully realized by deriving ex- 
pressions for the film Beyaold*s Htimber and mean heat transfer coef- 
ficient as follovs. 

Define a film Reynold *s Kumber, by 

» . ( -p- = ( vs)(^-yY) - ±L 

By makLing use of Equation (8) for P , ve nay vrite 
Rearranging the above expression, we have. 



( 19 ) 



( 20 ) 
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The original Musaelt relation gave cxaly i>oint values of the heat 
transfer coefficient and vas difficult to apply in an actual design 
study as it involved a function of the ccaidensate flow rate. This 
difficulty can be avoided, however, if a inean or average heat trans- 
fer coefficient is defined as, 






r»f 






Z At 



(25) 



By sijbstitutlng Eqixation (G) for P in Equation (2^) we have 



^Sk) 



Introducing the Prmidtl Kusber and rearranging. 



^ ^ _L fi 

T\2J f 2ot ' 7J\3j ^'[3 



2tx 



(i) 

(23) 

Hakir.g use of the dimensicailess terns 7 . , , and Lm , this reduces 

to 



k 1- ^ zr* 



I, 



j 2 ' 



(26) 



Equations (Ifi), (22), and (26) cccapletely define the results for 
the laminar region. For given values of the ditaennionleas sheer stz'ess, 
, either the mean heat transfer coefficient group, or the Kp^ 
can be plotted as a function of the dimensionless length of tlse plate, 

2 . This involves assming values of x^ end then solving for 2 from 
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E5^ia*iou (l8), frots (22), shkI frees £qusitioci 

(26). This has beea desie for various and for®s the the 

curves to the left of th® traasition points la Figures 7 to 12. 2he 
\iSQ of these curves in a design case vill he aiscaasad in a later eec 
tion, vhere the results of a turbulent fito as veil us the Iminar 
fila have beta derived, end are plotted in terms of these a«ie vari- 
ables. 
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3£CTim C - TBAmiTIOn 

For a liquid floving in film forai down a vertical surface it seeaso 
reasonable to postulate that there will be transition from laminar flow 
to turbulent flow at a critical Ileynolds number, in cajch the same man- 
ner that transition occurs in a pipe full of fluid. When we speah of a 
critical Reynolds nwsber it is to be roaessbered that this is not a def- 
inite, distinct point at which tiransition will occur in all cases. 
Bather it represents sobc sort of a mean value where the txarbixLcnce is 
sufficiently developed to saiterially affect the velocity profile within 
the f ila and where turbulent velocity fluctuatioas will be evident to 
the extent of affectizsg such characteristics of the flLm as the heat 
transfer coefficient. 

Before proceeding with a discussion of transition. It becomes ne- 
cessary to define a Reynolds number to give a cocason basis of coc^ieLri- 
eott. In the case under consideration it has become accepted to define 
the Reynolds amber in terms of the mass rate of flow per unit cross- 
sectional area, (w/O), the average fluid viscosity, ji, and an equiva- 
lent hydraulic diameter ccoaonly taicen as (^s/l), wi«sre Y is the unit 
width of the plate do^m which the fluid is flowing. In terms of these 
quantities, the Peynolds nmber may be defined as, 

/t 

Although there has been little experimental work ^Tith the specific 
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aim of dcterainiag the transition point of a freely falling vertical 
film, certain points may be inferred frcm heat transfer data. In ref- 
erences (12) and (15) on the basis of aaaevhat limited data it is rec- 
ocaaended that a Reynolds nmber of l800 be taken as the transition point 
for the case of zero vapor shear stress. It vill be shovn later that 
considerable doubt exists as to this val\^, however frota the present 
extent of knowledge no asore precise value can be predicted. Therefore, 
for the ifurposea of the analysis, this criterion for transiticxi will 
be adopted. One of the objectives of the experliaental work undertaken 
as part of this thesis will be to determine more adecpiate information 
on transition. 

Although we have adopted a critical Reynolds number of l500 for 
zero v&por shear stress, there remits the considerable problem of con- 
sidering how a shear stress may modify this value. Carpenter (5) has 
suggested that under the influence of a vapor shear stress, transition 
will occnar at a given value of wall stear stress. This laay be more 
readily vism,lized by considering the layer of fluid next to the wall. 

In laminar flow the sole effect of the outer layers cxi this inner layer 
will be a Bt»ar stress exerted on it. In effect, all the inner layer 
’’sees" of the outer layer is a transmitted shear stress. If we first 
consider the case of zero vapor flow, but increasing liiauid flow, this 
inner layer next to the wall will trip over to turbulent motion >4ien 
the Reynolds mrabcr attains the value l600. However, all the inner 
layer "knows" of this condition is a certain value of transmitted shear 



stress . 
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KoVjCon*id«rlng the case of a finite vapor shear stress, it is 
reasonable to assx»ae that the layer next to the wall vill chansp to 
turbulent flow whenever this "critical” shear stress is exerted on it, 
whether by outer layers of li^^uid or by aa outer laj/er of vapor. Tak- 
ing the wall shear stress, T© , as the loost convenient to which to ap- 
ply this criterion of ci'itical shear stress, wo may find its value frem 
the known case of transition at a Beynolds nuaber of iBOO with zero 
vapor shear stress. From Equatim (5)# T’ocan be expressed in terms 

of for w 0, 
o ' 



'Yo m -|- (p - p ) X 

go VO 

Similarly P , and hence Reynolds meber, 
in terms of x^ from Eqt»tion (8), 




may also be expressed 



(3r)t™> - leoo . V3 -^( 1 - ) xj 



Eliminating x^ between these two eqiaationn we are able to arrive at an 



expression for the critical wall shear stress, designated by ( 'T'o ). , 



( r©) 



tram 



J- 

% 



(23) 



For the more general case where ?£ 0, we apply our postuiaticai 
that transition vill still occur at this value of a critical wall shear 
stress derived above. From equation (5) we have an expression for Yo 
in the presence of a vapor shear stress. By Introdwing the vod.ue of 
wall shear stress obtained in Equation (23) and rearranging, the follow- 
ing expression for the dimensionless film thickness at transition can 
be obtained. 
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11.05 

Ti “ p.J7p~T 



1/5 




(£9) 



£iibstit\rting this value in Equation (22) give* an equation for the 
transiticsi FeyaolA* aiaa&ers in tersis of the <mly viarleble, . 



^ ^ ^traa ** ■’■ 0.667(1 ~ py/p) ( 

(50) 

This fuaetion is plotted in Figure 4 from & value of 0 to « 9*-» 
Inspection of Equation (J^) will reveal that whenever T»** > 

(i~p,yp)*'' ' 

* 

the dimensionless film tbichnees, x , viii be negative .*t t** nsition 

o 

Obviously this is lapcsslble physically and represents a ll-ilt to the 
applicability of tise prevlo\’,sly discussed criterion for traasition. 

In the higher ranges of vax^cr shear stress^ we isust, therefore, develop 
a mw critericn. 

For the hij^er values of it is quite possible that the film 
thickness at transition cq^isroacheB a limiting value, perhaps that of 
the ISBinar sublayer thickness which exists in ordinary pipe flow. 

Carpenter (3) suggests the following aagaltude for i^als thickness i 

• T ® 

Ordinarily a ssagnltude of 5 expresses the Isminer sublayer thiclmess . 
Actually an appropriate value it, soneyberu between 5 »ad per>-uk^ lo. 

The velus of Is arbitrarily chosen. 

If Equaticai (51) is solved for To and this expression substituted 
in Equation (5), the following cubic equation relating film thiciaiess 
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and vapor shear stress is obtaxrica. 



‘ 56 

J- - Py/P 






rJ' 



o t 



(^) 



For a given value of Ti- , may be found tram, this equation through 

« 

a trial and error process. With this value of and the given value 
of , the Keynolds nvaeber may be fomd through the use of Equation 

(22). The resulting plot is shovm in Figure U for values of If 9.2. 
It can be seen that the assimiption of transition s-t x^ » 6 is altaost 
equivalent to the assuaption of transiticn at a Reynolds miaber of about 
70 to 80 , there being little variation vf Peynolds number with vapor 
shear stress. 

This appears cn first ob&ervaticn to be a soiewbat startiagiy low 
Reynolds nucsber at which to expect transition, however, there is a log- 
ical basis for its selection. There is very little good data available 
to establish condltlcsis under which transition occurs when there is a 
significant vapor shear stress, A single point can be obtained from 
the enthanol data of Carpenter (5) by plotting his experimental points 
and noting the region of departure fraa the Uusselt laminar prediction. 
This ti'aasltion point at a ^ » 110 and Reyiiolds umber of 96.? agrees 
fairly veil with the transition Reynolds number of 71.3 predicted vmder 
the present analysis, m terms of x^"'’ the transiti<Mi data have the 
value of 6.9 which is close to tlie 6 of Equation (31). 

Although the criteria for transition ac portrayed by Figure ^ 
have not been adequately substantiated, the results are nevertheless 
quite reasonable euid are probably not serious Ij' In error. This entire 
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phenoaencm of trams It loa as affected by vapor shear stress Hfust be in- 
vestigated by obtaining nuch more experimental data. In spite of lack 
of adecjuate substantiation it cams reasormble to select the transition 

I 

criteria c.s shown in Figure h for xise in the present a^lysis. 
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SECTION D - TUI'BlILKIiT AKALISIS 



When condensate flow rates become large, it seems reasonable to 
predict that tiirbuleat flow will restilt, with considerable effect on 
the heat transfer coefficients. It was observed experimentally by 
Kirkbride (2) that, in the case of vapor condensing on the outside of 
vertical tubes, neasured heat transfer coefficients differed consider- 
ably frm those predicted by the Nusselt l«ainar analysis* In these 
investigations there was evidently no appreciable vapor shear stress. 

Colburn (8) reviewed the resvilts of Kirkbride and, by aasusing 
that the sole resistance to heat transfer lay in the Isainar sublayer, 
was able to derive heat transfer equations for a turbulent condensate 
layer. The Prandtl-Nikuradse velocity distribution was assvcved, and 
by selecting the thickness of the laminar sublayer at about « 10, 
a satisfactory correlation of the data was obtained. This equation, 
as presented by Colburn (8) is 




It is to be noted that the Prandtl nviaber enters here, vdjile it did 
not in the laminar analysis. This does not seem liuusual when one con- 
siders the correlation equations for heat transfer in pipes with no 
phase change. However, it should also be noted that this correlation 
VOS applied principally to Kirkbride *s data, with only small variation 
of Kp^ in the ra:ige of 2 to 5* 
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The ohove analysis ves extended by Carpenter and Colbux-n {5>9) 
to Include the effect of a vapor shear stress. In reference (3) heat 
transfer coefficients were investigated in a 5/0 inch copper tube 9 
feet long. Special effort was made to accurately raeasxire pressure drop, 
from which the magnitude of the vapor shear stress could be obtained. 

As the result of the expericieatal work and theoretical investigations, 
the following correlation equation for point valxies of heat transfer 
coefficients was obtained: 

= 0.0^3 [NPr) 

w 




By Integrating over the length of the tube, an equation for mean heat 
transfer coefficients was also obtained 



I , yX 

-= 0.0 io^ 

k 



(55) 



By the introduction of the diiiieasioaless vapor shear stress defined 
in Equation (17)# this equation may readily be converted into the form 

i. (N.y {r:f 

(56) 

The question of properly evaluating the vapci shear stress becomes 



ooraswhat complex when it is realized that, as condensation is taking 
place along the length of the tube, the vapor velocity is continually 
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decreasing. For convenience in Esaking calculations. Carpenter and Col- 
burn reccssaend taking a mean vapor shear stress, evaluated from the 
usviai plots of friction factor for one caaponent f lov at a vapor mass 



where is the value of G at the top of the tube, and Gg is its value 
at the bottom. It is felt that this represents scoewhat of an overs la- 
liflcatlon and that the friction factor should at least be evaluated 
frco the limited data cxi concurrent two phase flow available (11). In 
reference (ll) it is shown that there is considerable variation of 
friction factor from the snK>oth pipe data, mainly due to the fact that 
the liquid layer presents a ro^igh surface to the flowing vupoi*. 

One point especially worthy of conment in connection with the work 
of Caii»eater and Colburn is their consideration of the moaentum flux 
across the liquid -vapor interface. In their derivation they evaluate 
this tem as the product of the increment of condensate, ^ ^ and the 

main stream vapor velocity. This in effect assumes one dimensional 
flow of the vapor. In view of the feet that a vapor shear stress is 
also considered, it seams saae\rtiat contradictory to assisac, when eval- 
uating the mcmeiituia term, that a zero velocity gradient exists at the 
interface, and, wiien considering the vapor shear stress, to assme a 
finite velocity gradient. It seems more realistic to visualize the 



velocity. 
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layer of ve^or lnaaediately next to the interface as flowing with essen- 
tially the setae velocity as the interface itself. In this case it can 
be readily seen that such Eomcntuci terms will be aaall. It is throuch 
this line of reasoning that this nocentuo tern has been neglected in 
the present analysis. 

Inite recently Ssban (lO) has applied analogy calculations of the 
Prandtl-Karman type to a ttirbxilcnt condensate layer, assming the 
”miversal velocity distribution" of Rrandtl-^ikuradse . The improve- 
ment over previouis tuialyscs lies in the fact ttxat the resistance to 
heat transfer of the ccoiplete condensate layer Is considered, i-ather 
than dtist the resistance of the laainai’ sublayer. Seban*s analytical 
results extend over the considerable range of of 0.01 to 5 ajad 
agree well with the empirical results of Colburn (8) in the rang® 

Npj, »= 2 to 5* In the derivation Seben neglected the possible effects 
of a vapor shear stress. It is the purpose of the present analysis 
to extend £eban*s results to include this effect ai>d to complete the 
story by obtaining rusan heat transfer coefficients which include a 
laminar regicwi at the top of the plate, followed by a turbulent range 
over the lower regions. 

Turbulence has no effect on tl» force balance for a cojitrol vol- 
ume in the condensate film if ws think of the shear stress as an appar- 
ent rather than the Kewtonian viscous shear stress. Therefore, 

the resiilting force balance eqiiatlon will be identical with that for 
the laialnar case, and is given by Equations (2) and (3). The same lin- 
ear distribution of shear stress within the film will exist. 
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^ a -fi- (p - P^) (x^ - x) + (2) 

To “ “£"■ • Pv^ 

»o 

Siac« Kevtca*» law of vlacous shear will aot bold In the tm-tmlent 
region, it is Itaposslble to obtain analytically an expreasicn for the 
velocity distribution within the condensate fila. This analysis will 
be based on the assmption that the velocity distribution will be that 
of the Brcndtl-Sikuradsc formulation as given in refcrcfnce (12) . 
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The apparent shear stress, T » turbulent flow ray be expressed 
as the Bx«a of the viscoxis shear stress sad ti» eddy diffxisivlty tines 
the velocity gradient. 
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Conoldering the region in the film TJhere 50 < x^*", it may he 

eseujaed that • Solving for £kn > 



3. T- ^ 

= ~fT ^ 

Jx 



(40) 



Ibe term dv/dx may be evaluated for the region under consideration 
( JO <[ frota the asBumed universal velocity distribution. 



*» 5*5 + 2.5 In x^ 




The ccxitrol volvaac for the energy balance is shown in Figure 5 • 
Also shown is the tyi>e of temperature distribution to be expected for 
turbulent flow. 
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For the sliaded cofttrol rolurae it can be seen t^at there will be no ea- 
tnalpy cnanee for tne portion of condensate entering tnrougn tne upper 
surface and leaving through the lower surface. Further, it is apparent 
that the enthalpy change of the iacremeut of condensate will be essen- 
tially trAt frc8B saturated vapor to Bat\irated liquid. Therefore, an 
energy balance will be of the foraa. 



^ to . hfg 4 r 

As tine left side of the control surface is moved across the flla 
toward the wall, the above statements will hold true except for the 
thin layer next to the wall where the temperature distribution experi- 
ences a sharp change In slox>e. Dince this layer is thin, the effect of 
subcooling swy be ocsglected for the turb^llent film and q/A will equal 



(q/A)^, the heat transfer per unit area at the wall. Introducii.* 
this into the energy balai^ce, we have 



(JL ) . JL « h 4-C— (U2) 

' A 'o A fg dz 

The equation for heat trarnsfer across a turbulent film may be 
written in the form 

(Va I = (a + £,,) i± 

U 

where a = -h~ « thermal diffusion coefficient 

Nrr 

a eddy diffusion coefficient 

Through the application of the heat - monsentum transfer a^xalogy that 
<^1, » f end knowing the velocity gradients throughout the film 
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froia the assumed velocity distribution, v« are able to obtain an expres- 
sion for the t«nperat'ore gradients. By substituting for in 

Equation (^) and the velocity gradient from Equation (^1), ve obtain 

(Xo~^) 






fy 






T 



Substituting this expression in Equation (^5), the desired expression 
involving dt/dx and x os tl>s only variables Is obtained. 
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Solving for dt £Uid introducing, for convenience sake, the folboving 
notation, j? = =r ^ ^ ^ ^ (j~ A/^,) 
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Hhe condensate film .may be divided into the follovlrig regions: 
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figure 6 



Equation (kk) contmixis only tvo variable#, t atid x, einco we iaive sho-vci 
that q/k will be constant acroo» the file. This equation may then be 



integrated between the Units t to t 

^1 



^ and X, to X . After con- 
x_ i o 

o 



siderable algebraic saaiiipulatioa azid sixpllficaticai, the rc ,'Jlt of this 
integration across the film may be expressed as, 

2h 
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Equatlon (^5) gives &n express ioxi for the teaperature drop across the 
highly tiirhulent core of the condensate film. By utilizing the results 
of a previous analysis by V'on Karmaa (l4) for the temperature difference 
betveen the vail and the cuter edge of the buffer layer (x^ » 50), ve 
can obtain an expression for the total temperatiu’e difference froD the 
liquid -vapor Interface to the vail. Von Kartaan’s resxilts isay be put 
in the form 



(tx,~ f>Cp ]3(\-Hf)^- ^ 

iU) f ^ 



S t^Pf. f S Lh(lfS't^p^ 



(46) 

Adding Equations (45) and (46) to obtain the total tesaperattire drop 
across the film. 
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I^t the right hand side of Equation (47) equal Fg, then this ex- 
pression may be written 
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Por > 30 
o ^ 



p = + 3x^‘^ -h 2.3xp' In x^"") 



For 5 < <: 50 



r = (12.55 - 8.05 + 5 x^* la x^'') 



(>*9) 



Differentiating with respect to gives. 



For > 50 
o ^ 



=/-^(5.5 * 2.5 la x^*) 



4 Xj 



For 5 < C 30 



- /A' (-5*05 + 5 In x^"^) 
For convenience define as follows. 



For > 50 



3 5.5 + 2.5 In 



For 5 < V ^ 30 F, =• -5.05 + 5 In x ^ 
o 1 o 

Then Equation (50) may he written as 



For X > 50 
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JJl = ftJFi 



(50) 



(51) 



For 5 < X * < 50 4^ * 
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Equation 42 gives the results of the energy balance in the fora 

(-J9L) oe j r. 
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This may be written as 



(±\ - V II J)<>' 
\a L ' w H 



( 52 ) 



7 C- 















I 



O' < * *: ••'t 

I % ^ =- 

t>M»V • .* ••»- ’^•orA’- n 

I 

».«if:4 a* <l»#rf aOUAlf^iltitf 

' t ^ ’ I . 

^ > V ^ - 

I g AL 4 • «i*<-) ^ P i 

»P»*I I/I »« 

„.U .- •.< • X K •:> -^ 

'•.j/ {»t- ? ■« *f '>'.i^ < tis 

• # 

r.i • ■ (iVJ U 

=. liL 



*•* ^ 



.# > ^ 






^■*c Via '“ 






- 31 - 



Frcsa Equation ve have an expression for (q/a) as a ftmctioa 



of film thickness, vhich nay be r-e®jrranced in the form 



f _ At z’ Cf, (a:) 
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(55) 



Since it has been shown that (q/A) - (q/A)^, this term aay be elisslnated 



between Equations (52) and (55) • J^ubstitutlng for frcsa Equatim 



iL 

JXo^ 



(50a), we obtain an expression ccsitalning, as the only variables, , 



dx , and dz. 
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lliis may be integrated between the limits to x^ end z^ to Zg vhere 
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‘ot 



X' and z^ Indicate the values at the transition frcsa laain?-r to turbu- 



lent flow. Ills resulting equation may be put into the convenient form 




Equation (5^) uow gives the desired expression for the distance down 
the plate, z, as a function of flla thickness. In solving for this 

41 * 

distance, the value of x^ at the transition point must first be solved 
for frcci either Equation (29) or (32) and then the distance to the trena- 
itlon point, z^, found frcei Equation (l8). ITie distance fro© the 
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trajioitioa point to the point mider consideration may then be calculated 
from Equation (5i^) by evaluating the integral. This must be done nuraer- 
ically. 

Vhlle point values of the heat transfer coefficient as given in 
Equation (48) are of acad«aic Interest, they do not give aiiswers readily 
usable by the designer. To predict the heat transfer characteristics 
of a condenser or to enable rational design, it beccfflues necessary to 
predict mean heat transfer coefficients which can be applied to the 
condenser as a whole rather then at just one point. Accordingly It 
is possible to define this E«oa heat transfer coefficient in the semie 
manner as was done for the wholly laminar film. 

L = ^ =r ^ 

(t>, 

(55) 

It should be noted at this point that the substitution of q « w 
» r Y h - implicit in Equation (55) neglects the effect of subcooling 

Ig 

of the conderjjate film. As has been shown previously this is very near- 
ly true for the turbulent film, but it is not time for the laminar film. 
Since, in the present cose, we are considering both a laainar and a 
turbulent film, this will introduce setae error In the results. There 
will be no error In the lerainar curves provided an equivalent h^^ as 
defined in Equation (10) is substituted for '3^^* error involved 

in neglecting subcoollng may be somewhat evident in the region juet 
following troasltioa, but will becotac less important as more and more 
of the film is in tiirbulent motion. In most of the cases encountered 
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the lesainar film and transition region will extend only a few inches 
from the top of the plate and the neglecting of the eubcooiing in this 
small area should not Appreciably affect the overall results fo:^ plates 
whose lengths may approach scveiral feet. It is also felt that the un- 
certainties of actually determining the transition point will introduce 
errors of at least as great a magnitude as the subcooling error, and 
therefore this latter point should not be belabored. 

Since we have assumed a constant temperatujre difference, the de- 
nctainator of Equation (55) *sy be broken into two pajrts. 



The term (t - t ) z. applies to the laminar portion of the condensate 

X V V 



film, and the term (t - t )(L - z. ) applies to the turbulent portion. 

Xq w ^ 

Fraz rlquations (16) ai\d (l8), 



Equation (5^) gives the ejpression for the tvirbulent length. Pxxb- 
stituting these two values into Equation (56) yields the final equation 
for the mean heat transfer coefficient in terms of functions of the 
film thickness. 
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The Reynolds Kvraber at the bottota of the plate is of scsae interest 
as it is a direct measure of the aaiount of vapor that has been condensed. 

i/r 

This and taay be readily evaluated from Equations (^9) 

V 7 * 

»B. - -256 - X2 V * 10 xj lii x^* 

For 5 < ^ 

. 50.2 - 52.2 + 20 x^;; In 

Equations (5^^), (58), and (59) represent the results of the tur- 
bulent analysis in much the satae way that Equations (l8), (26), and (22) 
represent the results of the laainar aruslysii • The sesae fonas of plot- 
tiiag are therefore used, the tiurbulent plots merely being continuations 
of the leaninax curves. A description of the numerical evaluations of 
these equations is contained in the Appendix A. 
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SECTION E - DISCUGSIOa OF AliALyfICAL RESULTS 



The results of the aaolysis presented in this thesis are contained 
in Figures ^ through 12. These curves are baaed on a rational appi”oach, 
and. shew the effect of a vapor shear atreas cai heat transfer in laiainar 
and turbulent film condeasatim for a wide range of end Hp^. The 
use of these curves will be Illustrated later by a saaple problem. 

It is significant that the plotted values of meesn or average heat 
transfer coefficient Include the effect of a loialnar condensate film 
CO the upper part of the plate and a twbuleat film aa the lo>«r regions, 
if one exists. To use these plots it la not necessary to first inves- 
tigate whether or not turbialent flow exists, this deteralnaticn having 
already been made in the preparation of the plots. The values given 
are miversally applicable and are not dependent on tho geccsetry of 
the system or the flxiid properties except as they enter Into the deter- 
mination of the dimensionless quantities needed to enter the curves. 

Lnspection of Figures 7 througia 12 reva&l several interesting points 
concerning the raechanl«3 of film condensation. First, the significant 
effect of a vapor shear streos is apparent. This effect is more pro- 
novaaced in the lower range, giving an increase in heat transfer 
coefficients of several fold in sane cases. At higher the effect is 
less. The effect ca of the exact location of tiws transition point, 
as influenced by tlMK vapor shear stress, becomes less significant as 
the value of at the bottoai of the surface increases and is practically 




-I 



itt»ignificiuit for greater than 10,000. 

The effect of txirbulence in the f ila is illustrated by the change 
in slope following the transition points. It can be seen that the tur- 
bulence effect on the heat transfer coefficients is greater for increas- 
ing and for increasing Reynolds Suraber effect aight be 

expected, as larger will result in a greater degree of txirbxdence 

iiC 

and hence a greater portico of neat transfer by tins mechanise of tur- 
bulent lalxing. The Prandtl Ruaber effect is best seen by caaparlng 
the heat transfer plots for the very low Kp^ » 0.01, representative of 
liquid aetals, and the plot for the high « 10. The curves show that 
turbulence has a greater role in increasing the heat trarujfer character- 
istics for the hi^ier Wp„. In fact the low Prandtl maber cuirres for 
the turbiilent portion actually lie below esn extension of the laninar 
curves. Ihis may be explained once it io realised that for these liq- 
uids, the thermal conductivity is so large that it bcccaics the principal 
means of heat transfer, rather then turbulent nixing. Since the assuaed 
velocity distribution results in a thicker condensate film than would 
exist for laminar flow, it can be seen that this thicker film offers 
a greater resistance to heat flow through ccaiduction than would tlio 
laBslnar film. 

It would be desirable to cempexe the present analysis directly 
with cjqperlEiental results. However, experimental data on film conden- 
sation in which the vaqxDr shear stress has been accurately determined 
was not readily available. There does exist a correlation equation 
based on considerable experimental data, which has already been motioned. 



end ie given by Eq\iation (36). In order to caapare the anelytlcal 
results with this correlation eqmtloa, the curves in Figures 8, 10, 
end 12 have been cross plotted at a Heyaolds Staber of 5000 and ore 
presented as cxirves of mean heat transfer coefficient versus l/v . 

Theme plots are represented, by the solid lines in Figure 13 for various 
Prandtl Numbers. The correlation equation in the form given by Equa- 
tion (36) ^has also been plotted, and is sbcun by the dotted lines. As 
can be seen, the agreesnent is quite good for ■ 2 to 5> vhich lies 
vithin the range of 2 to 5 of the liquids actually tested and on vhich 
the correlation equation is based. There is considerable deviation out- 
side this range, as evidenced by the curves for * 10 and 0.01. This 
indicates that caution si^ruld be esployed vhen uaing the correlation 
e<^tion outside the range for '.ihich it vas proposed. 

Figure 13 also graphically illustrates the considerable effect 
that a vapor shear stress may have on the beat transfer coefficients. 

For exsKiple, the group for Np^ » 10 and ^ » 50 is about 1.00, or 
four tines the value of 0.25 given for » 0. 

In the force balance given by Equation (2) the mly pressure grad- 
ient in the z direction considered vas that due to the hei^t of the 
vapor and is represented by the A' tern. 'While this aasunption 
was valid for the case of a flat plate in an Infinite atmosphere vhich 
ve considered, it nay lose its validity vhen a vertical pipe is consid- 
ered. For pipe flov, the cxlstance of a vapor shear stress is intiaaately 
connected with a pressure gradient in the z direction. To include this 
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effect, it is nesessary to revrite the force balance. Equation (2) in 
the fom 

r- {.± f - ^) - r. -0 

( 60 ) 

To evaluate , cmsider a control voluae containing the vapor 

<y ^ 

in the core of a vertical pipe, as shown in Figure l4. 




FIGIBE Ih 

Considering the vapor flow an one dimensional and including the chaunge 
in ooaentua, the pressure gradient maj’- be evailuated from a force bal- 
ance on the control volme as. 
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J2. 3o Jz. 

( 61 ) 



where u represents the vapor velocity. 

2 

The weight rate of flow of the vapor, u n r^^^ , which my 



be differentiated to give, 

j ^ TT JjL 

Jz ^ Jz 

The weight rate of flow of the liquid condensate, w » 2 jc r^p i^ich 



may also be differentiated to give. 



Since 

<1 z 






- 2n IL 

iz Jiz 

, we may solve between these two equations for the 



velocity gradient, 

liL^ - A. 1£ Ji. 

JZ dz VV 

Substituting this eaqyressicaa in Equation (61) gives the pressure grad- 



ient as 



<lz 3o % n jz n y-c 

(62) 

Previously, in the case of the flat plate, we had taken the pres- 
sure gradient « (g/g^) Py* However, for the case of the vertical pipe 

condenser vs my take an equivalent vapor density 





JLAJUJI 

3 Kt Pi Jz 



a ia Ir. 

<3 Yi 



( 65 ) 



This equivalent vapor density, which takes into account the pressure 
Epradient for a vertical pipe, contains two terms which are dependent 
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on the radius of the condenaat* layer interface. For a given tijbe, 
the equations presented in the analysis can be evaluated in a manner 



the length of the tube, these must be determined from the given ccaxditions 
at each step in the stepvise maierlcal integration. It is apparent 
that the solutions so obtained vill apply to only the particular geom- 
etry for vhich they vere calculated and for the particular vapor vel- 
ocity which was considered. In this procedure, the advantage of the 
universality which was obtained for the flat plate will be lost. 

In Equation (65) the tenas containing the vj^r shear stress and 
the vapor velocity are subtractive and may quite possibly largely can- 
cel out. If the net magnitude of Equaticm (65) is email compared to 
the condensate density, p, it is then possible to apply the results 
obtained for the flat plate. 

It is worthwhile to Investigate the magnitvide of the pressvtre 
gradient correction for a specific case. Consider a 1 inch dlaraeter 
tube, in vhich mter is condensing. Satiirated vapor temperature will 
be taken as 120° F., and the tube wall temperature as 60° F. To take 
an extreme case, assume * 50. 

IJean film tsmperature, froa Equation (1), 



I p 

analogous to that for the flat plate. Since and u will vary down 



t^ » 120 - 5/^ (60) » 75° F 
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a point in the tube where x w 50« Then x « 500 end 

o ^ 



At 75° ?•» the Prondtl Kusber of water is 6.29. Proceeding in tlie Kjaa- 
uer outlined in the section on calculation methods, Fg may be evaluated* 



Fg * 57.6 

Then taking k of the condensate es .357 BTU/hr-ft •?F, 
h - 2660 BTU/hr -ft^ - °F. 

Sxjfcstituting h, *= 1CS6 BTU/lb, emd » 60^F., 

4^ - 2.59 Ib/hr - ft^ 

J2. 

Taking anooth pipe friction factors for illustrative jEHirpoees, 
f » .005 at a vepor Reynold’s Kvnber of 64,500. Since T/v< - t 
we may solve for tte vapor velocity, u, 
u » 1090 ft/sec 
("pe’v - - 64,500 

Substituting our computed valties in Equation (65), we are able to 
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evaluate , 
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4 2.34 - 21,82 » - 19.1^ Ib/ft^ 
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Although p^* is sbout l/5 of the condensate density, p, this Is 
not the loagnitude of the error involved by neglect. ing this tern. Re- 
ferring to Equation (bO), by continuing our pres«Kit cpcclfic exsmplc. 



The final crroi' which is iatix>duced in this specific case is of 
the order of 15.5^. These calculaticKis are included, not to exhaust- 
ively investigate the effect of the pressure gradient term, but rather 
to give an order of magnitisie of the eiTor which might be expected. 

The val\«53 chosen in the problem are felt to be representative except 
for the vapor velocity of 1090 ft/sec. This high value was deliberately 
picked to Indicate an vqjper limit to the error involved. In any prac- 
tical steam condenser, the vapor velocity will vsidoubtedly be much lower, 
with lower vapor sl»ar stresses. It may be concluded that in many cases, 
the eifect of the pressitrc gradient may be neglected, end the analyti- 
cal results for the flat plate laaj* be applied to vertical tube conden- 
sers with CCTifidence. 

In all of the calctxlations for Figures 7 through 12, the quantity 
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(I - p^p) t&k«a fui unity; ••g*, p^Jp Kaeill. It vae also assumed 

that acid it - t ) were unlfona along th« plate. The equatioiis 

o 

In the enalysis be x'eadlly evaluated accounting for any valite of 

(1 - p /p) and for any prescribed variations of (t - t ) or of Ta^ 

o 

along the length. The effects of these variations are fiHKrtioas of 
the particular goceietry and arc bard to generalise. 

OC the tsany assuaptions imder lying this analysis, most are sound 
and have considerable e:Q)eri*ental evidence to back up their validity. 
The two moot subject to q\Ksstion are the aasvraed velocity dlstributioix 

end the criteria for transition. It will be the purpose of the secood 

section of this thesis to Investigate the hydrodynaalcs of a falling 
vertical liquid fllsi under the influence of an appreciable sisear sti-ess 
at the interface. ISiese Investigations ere prelliainary esaly and will 
merely point out the direction to be followed in a more detailed exper- 
imental prograra. 




I 
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SSCTIOa F - SAMRE EgSIGI? FRQBIEMG 

Application of the results of the analysis Is conveniently 
made with the use of the ditaensioaleso plots of Figures 7 through 
12. Since aean or overall beat transfer coefficients are dealt 
with throughout, no point by point or stepwise calculations are 
necessary. The actual nethods involved in a design calculation 
are Illustrated by the following problems. 

The two problcsas selected are felt to be representative of 
those Ilyich night be mat in Indxistrial practice . The operating 
conditions and tube siaes chosen correspond to two actual cases 
where asmsxarenieats have been 'fcahen. The calculations ahcmld 
therefore serve the additional purpose of cccEparing the accuracy 
of the predict Iwis made by the present analysis. The first prob- 
lem utilizes the initial data for water nsi 191 of Carpeater’s 
experimental program, reference (3). The secend is based <xi an 
actual Industrial design as repcarted by Wuritcr, reference (15). 

The problasn is posed la the following aaaner. For a given 
condensate weight rate of flow per tube, it is desired to deter- 
mine the reqxiirod length of tube to condense the *smcrtmt speci- 
fied, with no uneonder.3ed vax.'or passing out of the bottcea of the 
tube. In addition, it la required to compute the mean condensing 
side heat transfer coefficients. 

In both cases the condensing vapor is steam at essentially 



atmospheric pressvire . The vail temperature has been token es 
the arithmetic average of the wall temx>er&ture8 I’eported in 
each of the references. 
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Ro. 1 


No. 2 




^be Chaa'acteristico : 


5/8”, 1^ BWG 


l’',17ByG 




Inside Diameter 


.459 


.884 


in. 


Cross-Sectional Area 


.00115 


.00426 


ft^ 


Tube wall temperature, t^ 


162®P 


182 .8® F 




Desired condensate flow rate, w 


50.1 


201 


Ib/hr 


Vapor mss velocity at top of tube,Gj^ ^5,600 


47,100 


Ib/hr-ft^ 


Since Gg^O, Eq. (57) gives 


25,500 


27,4oo 


Ib/hr-ft^ 


Specific voliate of saturated vapor 


at 14.7 psia, 


v» 26.8 ft^/lb 




Mean vapor velocity, 


188.2 


204 


ft/sec 


Condensate flow rate per unit periiaeter 






at bottom of tube, P »w/y 


4l8 


870 


Ib/hr-ft 


Average condensate flow rate Pt.. 


209 


455 


Ib/hr-ft 



Since water is the fluid being condensed, the ratio of surface 
tensions, s ^s = 1. 

The density of saturated liquid at 212® F., p * 59-7 

‘rhen 

>.o 



7.20 



5-5 

The viscosity of steam at 212'' P., ■ .0505 

f)C„ 

The Peynolds IJmber of the vapor, 29,700 66,600 

Utilizing the two caaponent flow friction factor data in reference 
(11), f ^ .0065 .0072 

rw f / 

f ^ — .1296 .17^ 



Ib/ft^ 

ft^/hr 

Ib/ft-hr 



Vapor shear stress 



lb/ft‘ 
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Mean Film Tempera tvire , free Eq. (l), 
Kinematic viecoeity at 



Ko.l Ko.2 

17 ^. >° 190 ® 

.599 X lO"'' .568 X 10“^ ft^/sec 
.790 X 10“* .7^1 X lO"^ ft 




Viscosity of cceidencate at t^, JH 

Ccaodensate Feynold* Kti^-jer at 

HV 

bottom of tube, N.-.^ » ^ 

Prandtl fKsRber of CcndenBcivte at t^, Np^ 

Enter Fig. 10 for N-, = 1 vith 

r; , get \ '■ 

Enter Fig. 12 for =• 10 vith 

ee. r:,get 

Interpolate to actual 

Thenaal Condwtivity of condensate 
at t^, k 

Mean heat transfer coefficient, h 

la 

t* - 

since t* • Cp / 3 I ... . 

Np*> hfjj w7 
solve for z 

Actual s, as determined by experiment, 
(references (5) and (I5) 



27.^ 


59-5 




,75^ X 10"^ 


X 10‘ 


Ib-scc/ft® 


1915 


4420 




2.26 


2.02 




.41 


.47 




.74 


.93 




.46 


.55 




.586 


.390 


BTU 

hr-ft-®F 


2250 


2470 




4l60 


8550 




3.6O 


10.4 


ft 


3.6 


11.7 


ft 



By ccs'siarifig calculated and actual tube lengths , it can be seen 
that the analytical results vrill predict the performance of a condenser 
quite closely. Probably the greatest source of error lies in the selec- 
tion of the com*ect average vapor sheer 



-C.!- 



■ fu P rr** , 



Si' 



5 c "^01 I *^fT. 



i f*' .'J^i I 

4i» 

A ^ 












'*e; 



rS**»^ t u . ^Qi > WT. 



'j*)i 

• i.«* 



^•fVi 



r^*i 






vrtt. 



«*“0 



4. -a 



OBi* 









"Xi) V-VMe 






& 



3 % >r 




« •« 



f J^f, • . 

"'‘(-1\ J 

KtJ^- 

i 4 « £a ' ^ 



‘.vr^^A 



dS*' 






ve./ 't« 



_ *y ^ 

^ *A ht\ 

» 

• V ^ J! 

<VT 

.*• f r j 






» 



« ftn* 



L .. '■''' -A ^2 ■ 

LVA-,) '>J ^,,1 V^n 

,j mtrm 4m V •• - ' f^ fci »4 

i) b t 



s eJ ii' f 



• * r , • V. »-w •.«. ^ 



• m 4 t %4 

“J - xJ^mrit • WX^ 



M« %» .«Mf 



- 54 - 



str«S8 ciid tube wall t«aap«re.ttire . If it vas deeired to predict 
perfonaance within rery elo«e tolerances, the niEierical inte- 
gration procedure described in AppezKlix A could be reworked witij 
the known variations of ^ t and included in a stepwise appirox- 
imatiem. However, this vouM involve a rather tedious and lengthy 
set of calculations and, for Bc«t eases, it is felt that the cm- 
stant averace values as used above will provide resui^ of s\jf- 
ficient accuracy. 
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part II 



EXPE^IMEIjTAL LWESTIGATION OF CO«CUHHEI?T 
LIQUID-VAPCgt n.OW IN A V3RTICAL TUBE 
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PAj’^T II 



BSCTION A - LriERArJgS SURVEY 



Part I of this the*!* is coticemed vlth tlie developraeat of ana- 
lytical expressions predicting heat transfer ccjefficients of a 
falling fllxa condenser. Use has been made of the heat Eomc-ntm trans- 
fer onalogy, so that the validity of the predicted heat trensfer coef- 
ficients is dependent upon a knowledge of the hydrodynsnlcs of a falling 
filjB. Assun^tions have been aode rej^rding the behavior of such a film, 
based on the best knowledge available today. It is the purpose of 
this section, therefore, to give a brief review of tiie literature cov- 
ering the subject of falling lltjuid films. 

Honrad and Badger (15) discuss omet of turbulence in a falling 
liquid film. Their ccaicliwion, based oa the work of other iavestigators, 
is that the critical Ecynolds ffuBber is in the range of l800. They 
also cite data of Jakob and Erk (19) as showing that transition say 
occur at lower N„ in the presence of high values of velocity. 

In reference (17), Rebbard and Badger neasured film coefficients 
of heat transfer on a 1 inch O.D. vertical tuba condenser 12 feet long. 
Beynoldo Kmbers of the film covered the range 500 to 5000. Results 
are coaslstently hi^jer than thc«e predicted by the Eusselt lorainar 
equations, and generally lie in a line parallel to the predicted line. 

Ro aseasureaaente of the vapor sheer stress were made, and it is felt 
by this author that the dlecrepeacies quite pcoalbly be explained 
by this effect, Ko information as to the onset of turbulcait flow can 
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be deduced frcsn these results. 

Kirlfcride, reference (2), presents heat transfer data of several 
invest iGstors on liqviid films flowing down vert.ical walls. Tliis data 
covers the Reyaolda Rmber range of 700 to 20,000. By the aarked dif- 
ference In slope of the correlatin^'j line fro« the Russelt laminar pre- 
diction, turbulent flov of the film is definitely indicated. Ho>«ver, 
there is little data for the lover 11^^, so it is not possible to direct- 
ly infer the location of the transition point. The best correlating 
line through the data Intersccta the theoretical Husselt prediction 
at a Hp^ of about 1100. The vapor shear stress was not considered. 

In addition to considering heat transfer coefficients, Kirkbride 
(2) presents data esn seasured flla thicknesses for less then 2000. 
Excellent agreeaent with predicted thicknesses is obtained for the lower 
Hp^, hovflever, deviations are observed for the higher Surface waves 

are given as the explanation for the discrepancies. 

Cooper, Drew, and McAdaais (l8) present isotaenaal data of six 
observers dealing with ll<3;uidB flowing down the wetted smooth siu’face 
of vertical towers or over flat plates at slight inclinations. Film 
thicknesses were aeasured either by a microQctsr, visually observing 
contact with the film, or by suddenly stopping the flow and measuring 
the qimtity that drains out of the test section. In a plot of the 
Fanning frictiesa factor versus scatter of the data frem the lami- 
nar equation starts at e. of about 1000, with the best correlating 
curve indicating transition at about 2100. Bata were not available 
for dsteraining the possible effects of gas velocity on the liquid flow. 













r 



•w 



I 



j 






jptfUM 



-5a- 



Considera2:ile expericeiital vork has been done on non -condensing 
heat exchongers ©sploying vertical tubes, in vhich the fluid flovs in 
film form dovn the tube vails, references (20) to (26). Bays in ref- 
erence (22) has plotted heat transfer coefficients versus A break 

in the curve at a of 1300 or ihoo raay possibly indicate transition 
of the film from laminar to turbtilent flow, however. Bays recoamends 
a value of 1000. Leavitt’s data (25) for water extends from a of 
1000 up, with no indication of transition as low as 1000, althouglx 
considerable scatter of the data is observed below !f„ » 1700. The 

data of Fentress (2k) covers the ReyiKslds Ehsaber range of 85 O to 20,000, 
with no break observed vdxich might indicate transition. In ail of these 
references, no vapor shear stress was present and no attempts were isade 
to measure film thicknesses. 

Dukler In reference (2?) has used a film thickness measuring tech- 
nique employing a capacitance measurement. This device determines film 
thickness by the difference in dielectric eexustants of air and water 
between two plates of a capacitor. One of the advantages of the tech- 
nique is that it averages out surface waves and irregularities. All 
measurements were made on the film flowing down a vertical, flat plate. 
Film thickness readings agree well with the l^xisselt prediction in the 
lower range and with the prediction from the "universal velocity 
distribution" In the higher range. As tlie technique is somewhat 
difficult to apply in a round tvibe, no investigations of a vapor shear 
stress effect wro mdertaken . For the flat plate with no vapor shear 
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streeSf transition is iridieated at a of 1000. 

It is apparent that considerable variation exists in the results 
of the many investigators cited here. On the basis of the vorks cov*' 
ex'ed in this survey, no dif inite conclusions can be reached as to the 
onset of twbulonce in a liquid f ila flowing downward under the influ- 
ence of gravity. Furtbemore; m> Investigations have been »ade to de- 
termine the effects an appreciable vapor shear stress. More isnovX- 
edge in this directicxi is urgently needed. This, therefore, vill be 
tlie directioa to be taken in the expertsMsntal pert of the thesis. 



f 
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SITTIOK B » £XP1;HIMST!!TAL £QUIR<ERT 

’Thus purpcwie of the expericsental program io to inveatigjatc 
hydrodynaeiics of a liquid filR flovlns dovn\«rd along a vertical eur- 
foce mder the Inf XiMsoce of a shear stress exerted on the cuter face 
of the liquid. In order to obtaizx quantitative 2 *« 6 ults It la neces- 
sorj' to be able to directly oeasure the following quantities; 

(1) Ihe liquid flow rata 

(2) The liquid taeperat’ore 

(5) The ffiagnitude of the interfacial shear stresa 

(^) The liquid fila thichaess 

(5) In ladditica, it would be desirable to be able to obtain 
oeasurementa \&ich would give quantitative Infonaation as 
to friction factors for two coK^ment flow 

r 

Althowigh a two dlaaaaioual set-^, with the liquid flowing down 
a flat pXate^ was considered, the necessity of accurately dcterainlng 
the Interfacial shear stress lad to the choice of a vertical tube for 
the test section. By introducing the liquid along the tube vail, it 
flows down the walla in film form. A flow of oii* in the core of the 
tiibe stoulntes the vapor flow in a vertical tube condenser and supplies 
the ei»ar stress on the outer face the liquid. By assuming no 
pressure drop across the liquid fi3at, this shear stress can be easily 
measured by presswe drop readings along the length of the tube. 

The test section is shown in Figure 15 and is a one inch diaaeter 
tube cppi*oxl«ately 5 feet long. Preesixre taps are installed at one 
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foot iatervsLis along its length. In oi*der to prevent arcii^ over of 
the liquid at high flow rate©, a «aooth convergent nozale is placed 
at the top of the tube. Surrounding the nozzle is a oettllng tank 
to daop out any turbulence or fluctuatloas which aight exist in the 
liquid as it is discharged from the st^ly systea. fhe nozzle acts 
as a veir in the settling tank to evenly distribute the liquid along 
the vNills of the tube. The entire test sectim, nozzle, and settling 
tazik are constructed of luclte to permit visual ebservatlon of flow 
cliaracter IstlCG . 

Six probe ports are installed in the test section. These are 
arranged in two groups of three each, spaced evenly arotsnd the cir- 
cmference of the tube. Details of the probe ports are shown in 
Figure l6. Film thickness readings are made with the sundifled depth 
gage micrometer, shown la place. i4etal probe plugs, whose inner ends 
are fliwjh with the test section wall, are installed to provide elec- 
trical continuity from tljo alcrosseter probe through the liquid film. 

The microBieter probe is a standard depth gage nlcrocieter with a 



except for the extreme tip, are insulated. This probe is shown in 



Fifjure 17. To obtain film thickness I’eadings, the probe is inserted 



blocks provided. A aero reading is obtained by advancing the micro- 
meter until it touches the far wall of the test section. The pi-obe 



is retracted end a sensitive oiaaseter connected between the end of 
the nicroaoter and csie of the adjacent port plugs. The microtaeter 




shaft. The shaft and the extension, 

I * 



in thfl port, with the shoaldcis of the micresjaeter held fimly on the 
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ic then •dvuaced until electz'icftl continuity i» cstebllshed by the 
probe tip touching the outer edge of the liquid film. Ihla contin- 
uity 1» easily detected by the ohaaeter. The insulation ao. the shaft 
of the probe prevents completion of the circuit by the film on the 
nesu: vail of the tube^ through vhich the probe extmads. 

Pressure readings in the test section were made with a standard 
U tiihe nenooieter vlth air in the tubing between the rxaooneter end the 
pressure t^. Surface tension forces at the interface between the 
liquid on the vails of the test section and the air in the tublTtg, 
plus the teadancy of the liquid to be drawn froG the test section in- 
to the tubing wAzen changing the presexurc level, stede tiieae readiiigs 
scmewhat difficult, yncertalnties of measured pressure levels were 
(mly evident at low air flow rates, vtoen the surface tension forces 
were of a ccc&parable magnitude to the pressure foirces. Modifications 
have been made to correct this by using a single manoseter tube, with 
the connecting tubing completely filled with liquid. This requlreo 
that a z«ro reading be tshcn with no flow in the test section, but 
eliminates all interfaces in the pressure measuring appeuratus. How- 
ever, ail data presented here were zaade with the original system. 
PresBixre readings plotted as straight lines for all but the very low 
air flow rates, and ajqpear to be fairly reliable. 

The test section discharges to a closed separating tshk, whose 
purpose is to separate the liquid from the air. The air suction line 
is connected to the top of the tank through a standard pipe Tee inside 
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the tarJc. The air flov oust traverse several s.^arp ciar.ges direction, 
thereby separatiiig a..y liquid droplets which m.%ht be entrained. A 
liquid recirculating pump is co;irected to ttie bottos of tiie separating 
tauK. through a flexible nose. 

Jio atteept was made to ueasurc tne air flow rat>r* in tJiese prelim- 
inary investigations. Gucn laeasiirements could be easily made by instal- 
ling a rotaaseter or seme positive displacetie it measuring aevice in tne 
air suction line. It is felt tivat tcese measuresaents, wriiich would en- 
able detera^atiou of friction factors, wuld be highly desirable for 
future experimental work. 

The liquid flow rate is neasured by an orifice iiastalled In tie 
supply line. The pressure drop across this orifice is measured by a 
water over mercury laaraaeter. This system was calibrated In place by 
time - welgx>t saeasurenents . Calibration data and curves are given in 
Appendix B. Liquid temperature readings are obtained from a tbemcmeter 
placed in the settling tank at the top of the test section. 

Auxiliary equipment consists of a pvKp to recirculate the liquid 
through the system, and necessary piping. For the preliminary experi- 
menteil program included in this thesis, water was the only liquid tested 
and therefore, connections were made to a water main as a source. A 
line sketch of the equipment is shown in Figure 10. Figure I 9 is a 
photograph of the set-i^, shoving the test section assanbly and sepa- 
rating tank on the right, with the micrcoeter pi’obe in place. 
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Slnce it %»as ar.tlcipated that vibrations would quite possibly be 
critical in affecting the steadiness of the flow, flexible hoses were 
installed in the pipii^ from the separating tank to the puap and in 
the supply piping Just before the settling taiUt. Ihe test section was 
ttaxA isolated from vibrations in the piping system and pump. 
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SECTIOH C - FIU4 THICK?n&S3 



Film thickness tteasureseats have been made la the asnner described 
in Secticai B, Part II. Since ve are dealing with a round tube rather 
a flat plate, it woixld be worth\diile to investigate the equations 
applicable to this case and conpare then with those derived for the 
flat plate la earlier eectims* 

Consider the control voluae, idiose cross section is shown shaded 
in Figure 20. This control volume is an annular ring around the tiabe. 





Since there is no condensation, SKxaentvea terete nay be neglected. 



and the force balance for the control volirae is given as 
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T^TTkc/z -X^^TTn H =0 

( 64 ) 

After S0G!G siiu.pliflcs,tion ccsid rcenTangasaenti this reduces to 



t' - X^T “T“ -J il y- 



< 65 ) 



Since mooentuii terms do not enter, it is possible to use the ex- 
pression for of Eqixation (61) with the second term * 0 



Substituting this in Eqiyition (o5), we get 



( 66 ) 



r- 1 (/’-^~) *- F 

( 67 ) 

la solving for dp/ds fress conditions in the vapor core of the tube 
and then applying this dp/dz to the esjuations for shear stress in the 
liquid film, tlie asswjption has been mde that no change in the pres- 
sure gradient occurs in a rfl.*lal direction. This assisaptioa Is quite 
good, and aOmost exact for thin liquid films. 

m order to cempare the distribution of sheai- stress given by 
Equation (6?) with that pz-eviotisly derived for the flat plate, it is 
necessary to substitute ead r * - x. These substitu- 

tions ore Identically true from the geometry and ®ez*ely serve to shift 
the zero axis froro the tube centerline to the tube wall. Making the 
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siibstituticsis and sls^lifylixg, ve have 







^ Vp-x 
y«~Xo 



( 68 ) 



If tt» liquid filia is thin, r » x and r » x . For this condi- 

' o o"^ o 

tlcaa conaiderohle simplificatlcai nay be made, resulting In the follov- 
ing shear stress distribution. 



This expression is identical to Equation (2), which was derived for 
the flat plate, and forms the basis for the development of tlte predicted 
heat transfer coefficients. 

In the experljnental program, ths ffiaxisum film thickness measured 

was of the order of .025 inches. This represents of the tube radius. 

3y sxjbstltutina: these values in Bquaticxi (65), it is easily shown that 

the nsaximmn error which is introdiiced by the simplifying assisaptioa 

that the film is thin is only 5.5%« This maxlEm error will occvtr when 

considering the wall shear stress and will only be apparent when the 

vapor shear stress forces are large compared to the gravity forces. 

In this case, the film thicknesses are much analler than .025 Ixxches. 

Having verified that tije flat plate equations foi- shear stress 

are approximately true for the case of a vertical tube, it is now pos~ 

Bible to utilize the flat plate equations to predict the film thick- 

ness as a functlcai of . For laminar flow, Equation (22) applies. 

♦ 

For a given Feynolds number, by assuming various values of x^ and 



T= x (/-A)(x.-x) i- 



(69) 
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r / + * 

„ V , we are able to plot as a fmction of . 

Similarly, for the case of tvtrbulwit flow in the liquid film, for a 
given KeynoMB Kuaber, a trial and error ooluticxi of Equation (49) 
vill yield a value of . From the definition of "quation (j3), 
the wall shear stress can be found in terms of the film thickness, x . 
By equating this expressicai for wall shear stress to that from the 
force balance. Equation (5), a forraula for x^ in terms of Is 

obtained. This iias been done for both lasainsr and turbulent films at 
the Reynolds Bumbars for which fito thickness readings were obtained, 
and forms the basis for the theoretical curves shown on the plots of 
measured data. 

Tlie procedure used in obtaining meastured values of film thickness 
is as follows. Liquid was recirculated witii tiie pxsap, csnd the flow 
rate adjusted to give the Feynolds IJUtaber desired. The purap w&a vised 
all of these runs, rather than water directly from the water main 
as it was found that fluctuations in city \mter pressure made it .impos- 
sible to obtain steady flow x-ates. With the water flow rate constant, 
the throttling valve on the air suction line was adjusted to give var- 
ious levels of presoxxre dre^ through the tost section. At each of 
these valxiBS of pressure drop, film thickness readings were taken with 
the micrometer probe, following the procedure described in Section B, 
Part II. Before each series of i*uns, film thickness readings were 
taken through all three of the probe ports in the grox;^ being worked 
with to insure that film thieloiesses wexe equal around the circurfereace 
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of the tii>e and that the flow vas evenljr distributed. The raw data 
and ccmputations for these njns are included in Appendix C and ore 
shown plotted in Figux'es 21 through 27. 

The static pressure level at points in the test section was plot- 
ted against distance, and a straight line curve drawn through these 
points. The slope of this curve vos then taken as the presaiu*e grad- 
ient in the z -direct ion. This pressure gradient could tixm be easily 
converted to the vapor si«ar stress, , frees the formula 




( 70 ) 



where is the raditis of the tube 
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DISClB8I0fi OF FID4 THICKJIJKS MEASURISStlTS 



Inspection of t*iS plotted results of the film thictoess measxire- 

mento. Figures 21 to 27, reveals several interesting facte. First, 

* ^ ^ 

the shapes of the curves of versus Lm correspond closely vith 
the Bhaj>eB of the theoretical curves. In general, the values of meas- 
ured film thicimess lie below the predicted values, this deviation be- 
ccEilag more noticeable at higtier liquid Reynolds Numbers. The neasxjred 
values for L/.- „ q are consistently higher then pi'edicted. At high 

>y*‘ * 

values of , experimental values of agree well vith theoreti- 

cal values. There appears to be rather sharp change in curvature of 
the measxircd values at values of in the neighborhood of 5 to 10. 

For t- a 0, visual observations revealed tlxit rather largo sur- 
face waves were present on the interface of tiie liquid layer. These 
waves made film thicloiess measurements, using the micremeter probe and 
electrical continuity technique, semewhat uncertain. The otaaeter 
registered alternate breafeing cad making of contact between these sur- 
face waves and the probe. In taking the readings, an attompt was made 
to obtain a meon wave heigiit reading, but this procedure is admittedly 
a crude one, as surface t«isicn forces tended to maintain contact be- 
tween the probe tip and tlae liquid surface. Therefore, it is not sur- 
prising that these readings were high, it being believed that thick- 
nesses, as Eeasured, were actually somewhat nearer the crests of the 
waves than the mean height. 

For all val\jes of greater than about 5, the large, slow moving 
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Btirfaco vav«o were replaced by arueb smaller and more rapidly movins 
irregularities. Indications on the ohpjmeter gave quite sharp cut-off 
points for contact of the probe with the surface. It is felt that film 
thickness readings for the higher values of are accxirate to within 
.(3005 inches. This is l>ased on the sharp cut-off points and the con- 
sistency of the readings chi several rune under identical condltiosis . 

* 

The fact that msastired values of lie below the predicted cur- 
ves leads to the conclusion that the velocity distribution is actmlly 
sccuething other than that aesuaed. Although the theoreticad. turbulent 
curves lie above the leainar curves, it is believed that, for the higher 
values of , the film is essentially in turbulent notion. This con- 
clusion is based on the eppearonce of the film and on sooe preliminary 
results obtained with a hot wire onenonieter in the liquid film. The 
fact that the velocity distribution may be something other tlian the 
ruiiversal distribution of Prandtl-Jftkuradse is not altogether srirprls- 
ing. The use of this velocity distribution is certainly en over exten- 
sion of their work, ^ich was derived frcaa measurements of full pipe 
flow. However, the actual values of film thickness are at least with- 
in a few percent of the predicted results. 

It is believed that the change in curvature of the measured values 
of fil 2 i thickness is due to the damping out of the large surface mves. 
This corresponds to visual obscrvatloas . The general character of the 
flow, as observed through tlie transparent test section was as follows; 
For « Of the flow was mainly ciwir with relatively lax-ge and reg- 

ular svcrfajce waves. The frequency of these waves was low and fairly 
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constaat, being of the order of 8 cycles per second. As the air flow 
was increased, /these largs surface waves were dataped out until, at 
around Ta^ « 5# the flou appeared quite clear and transparent. For 
yet higher vjiilues of , verj' simll surface waves were formed, with 

a x'clatively high fx*equoncy. It is estlasated that the height of these 
waves was less 'than .0005 inches. It is felt that these asaall waves 
are directly connected with the onset of turbulence in the liquid film. 
It is to be noted in this respect, tliat tlujse aaall waves could be 
attained without any air flow at all, but required very high liquid 
Reynolds Syfabers. No measurcBtents of this point were obtained, the 
flow rate being beyond the capenity of the liquid flow sxancnxeter. It 
is estlnal|ted that liquid Reynolds Ifusibers of the order of nacpitude 
of 6,000 to 8,000 are required to attalxa this condition. 

Cto the basis of the film thickness aeasurementB , no definite ia- 

f 

fortaatlcxi could be ded^^ied as to the onset of turbulence, as affected 
by a vaqjor shear stress. However, the equatims predicting film thick- 
ness pro\'«d to be approxljnately correct. It may be concluded from, this 
infornjation that the assumed velocity distrlbutlcai is adequate for the 
purpose of predicting heat transfer coefficients. It may also be cm- 
cluded that the equations derived for a flat plate will be equally 
applicable to vertical tubes, as long as the film thlclciess is asaall 
ccnq>ared to tlxe tube radius. 
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SECTIQK D - WSimimnY BET.UIJ.’S (F IIOTT w&E AKEMCMJTEr! 



In Paurt I t>#c criteria fcr datura Iniag the traiisitim point cC 
the liquid fila from lawinar to tiirbulent flow were propc«Bed. These 
criteria, althooi^ based on some experimental work, are far fx^cm beiixg 
verified, especially in. the region of ^preciable vapor s^xear stresses. 
In fact, no experimental work has been found that con shed light on 
this phenoBMsnon of trarjiition, as affected by a vapor shear stress. 

This section will discuss the prelimiaary work for deteraining the tran- 
Blticffi point by means of & hot wire aaiaacoseter . No quaxitltative Infor- 
mation has been efctained, but work along these lines is currently in 
progress and should becocie available in the near future, reference (l6). 

The instxnaaent is a standard hot wire aaesaosaeter, manufactured 
by the Flow Corporation of Cambridge, Jtoasachusetts . It is being used 
In cenpressor blade cascade studies in the Gas Turbine laboratory at 
M. I. T. The elesftent is a .0005 inch diaeioter tungsten wire, supported 
by a long shaft. By inserting the shaft through one of the probe ports 
on the test secticax, the wire can be introduced into the llcjuld film 
on the far side of the test section. With currents of approxlta&tcly 
150 milllsaspa, good indications can be obtained eax an oscilloscope. With 
the relatively high cnarrent, the wire burns cut almost Instantly if ex- 
posed to air. This has been the source of scaie difficulty, end it laay 
be necessoary to modify the wire element or lower the current to pre- 
vent this. 

At low liqsjiid Keyaolda f?ijEbere, tlic flow appears to be quite 
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siysceptible to unsteadiness caused by vibrations. PreUfalnai-y runs 
vcre aade vhile a large air coraiiressor was runnlns in the same build- 
ing, ai>proximately 100 feet auey. Very regular velocity fluctuations 
were observed at lov Bejuclds Ifutsbers, which died out as the liquid 
flow was increased. The runs were then repeated at night, when the 
coEcpresBor was not running, and these fluctuations were absent. 

For the case of no shear stx’ess, laa in a r flow in the liquid film 
existed up to Peynolds !?UBb€r of 6,000 to 8,000. A trip wire was in- 
serted around the noszle throat at the entrance to the test section. 
With this trip wire, occasional turbulent fluctuations were observed 
at tlw* flow rate tlmt the liquid flow riaanoBeter could meiisure. 

It is believed that stable laminar flow exists up to Kcyaolds fiteabers 
of 6,000. This observation has at present not been verified. 

With air flowing la the test secticsi, a consistently turbiaent 
flow couM be attained with lower liquid Reynolds Sumbers. However, 
as was evident from the film thickness measurements, the film is ex- 
trasiely thin for the regions of appreciable vapor shear stress, and 
difficulty was encountered In Keeping the hot wire completely in the 
vmter film. It say be necessary to use a more viscous fluid in order 
to get film thicknesses of sijfficient mag^iitude with which to work. 

On the basis of the very sketchy work done with this technique, 
little of a quantitative ixature can be eonclxxled. However, tJie tech- 
nique appears to be a powerful ojie. Turbixlent indications arat quite 
marked on the oscilloscope trace, although a lower bscKgromid noise 
level would be desireable to enable ao3i*e exact determination of tlwi 
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traaaiticHi fr«B ccsapletely Inalnar cooditione. It is felt that, once 
the minor problems have been solved, the hot vjire aaeaoBctci' v/ill yield 
the answer to the transition questicai. 
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SgCTIOS E - FOP FUTIEX EXPHiH^nkL 



In view of the present Iscic of knowledge conceixiing the behavior 
of a falling liquid film from a fluid dynacics point of view, it is 
reccBEaended that the fiirst experhaentai efforts be directed along this 
chennel. 

Specifically, it is felt that it would be desirable to fully in- 
vestigate the transition criteria. This should be dene with at least 
three liquids of widely varying viscosity. Consideraticai should be 
given to varying surface tension also. This receaxaeadatioa is based 
on the fckct tliat, in this type of fluid flow, these forces may quite 
possibly be the only significant eaxes in the r*adial directicxi. It is 
worthy of coiament that a siirface tension ratio was the correlating 
factor in the fricticai factors preswited in reference (11). 

In workiiig with other fluids, it would be desirable to verify 
the fil« thickness meastirements presented in this thesis. By gather- 
ing sufficient infoitaation, equaticais taore adequately describing the 
velocity distribution in the turbulent regions should be developed. 

Much more information concering friction factors in two phene 
cocurrent flow is needed. As previously pointed out, this can be 
easily obtained from the present set-v?? by adding a aster to measure 
air flow rates. 

After the ijydrodynaRiics of the liquid film have been sufficiently 
defined, the analysis should be modified to incltide any new or better 
infoimation that isny nave been discovered. Finally the analysis should 
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be directly chocked by Investli^atlon of &a expcriaental condenser, in 
\diich provisions should be laade to deterainc point values and raeon 
values of the condensing side heat transfer coefficients. 
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APPE9DIX A - mMSB.TCfiL EVALUATIOH OP AflALYTIC £'4llATIO;€ 



Laalnar Portion 

For this portion of the condensate film, the Bignificant Eqiiations 
are (l3), (22), snd (26). 



2* = (/(') 3 T-" (A.*J 



(IB) 



jjs . .... = i (x:) 



ar:u*f 



-+• 



1 =^ ^ X rf Ud 

/i U/ ^ z* 



( 22 ) 



( 26 ) 



fy’* * 

For given values of L'- assume various values of x lower than the 

o 



transition point, which is given by either Equation (29) or (52). For 

♦ * 1 
, both 2 and 

k (itf’ 



the assumed values of , both 2 and evaluated 



directly. 



may now be found frota the given values and 



the ccaputed values of z . 

Turbulent Portion 

For coeputation purposes. Equation (53) aay be rearranged into 
the more ccxivenient foita 



ib fiL"f - N,,. 

/f bv ■ (2*-zt‘; 

(63a) 



The evaluation of this equation entails & nuaericai int^ration pro- 
cedure and is best followed by listing the steps involved. 
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step 

Step 

Step 



Step 

Step 

Step 

Step 

Step 



Step 

Step 

Step 

Step 



1 - For given values of aad , etssumc varicue values 

FT 

* 

of gi'eater than that at the transition p/oint. 

2 - Evaluate M frean Eijuaticsi (^5a). 

5 - Evaluate *■ )(l - p^yp) . At this 

point in the ccsiputaticxis it vas necessary to assure that 



(l - P^P) » 1* 

^ - Calciilate Fg from Equation (^7). 

5 “ Calculate froo Equation (51). 

6 “ Form the product, F^Fg 

7 - Plot this p 2 *oduct aa a function of 

3 - At various values of , take the area isider this curve 
with a planiaeter or by other suitable graphical Bieans. 



This area is taken from the transition point, to 



the selected values of x 



. * * ii 

9 - Calculate (s - z* ) » 










10 - Add to the result of step 9 the distance to the trans 

■» 

ition point, , foxusd in the laminar calculations. 

11 - Calculate from Eq;uation (59)* 

12 - Nov ha. = Nji^ 

k l3/ 7* 



* 

Where z is the result of ctep 10. 

Step 13 - Plot versus a* end There will result a 

double series of curves with Tj* and as pareneters. 

Ft 

The results of the above calculatiaas ai-e contained in Table I and are 



shown plotted in Pigures 7 through 12. 
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APPEKDIX C - DATA AND CALCULATIONS FOR Fllil THICKNESS 



Theoretical Predictia^^o of Film Thickness 

laalnar - The laainar predictions result directly from Equation 
(22), which nay be rearranged in the fom 

'V'* _ Nfie. " y (^») 

x(x.T 



For a given the plot is made by asstsalt^g values of a»-id solving 

tor L/y- 

Turb u lent - The turbulent predictions result froK tiie followltig 

fee 



procedure; First, consider Equation (**9), which relates and 



For > 50 



^ » Np^/4 » (-64 ^ 5 + 2.5 x^" In x^^^) 

For 5 < x^^ < 50 



X - Np^/4 . (12.55 - 8.05 x/ . 5 x/ in 



For a given T!p^, this equation is solved by trial and error for . 
Now, by eliminating To between the ‘'cfi..*ition of x^’, Equation (38), 
and the force balance. Equation (3), Tv vmy be expressed in tenas 

-f * 

of X and , 
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With the value of x^^ determined by the Np^, and by assuming various 

valties of X , can. be solved for. 
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The above pi*ocedxires for both the laainar and turbulent films vlll 

yield the predicted film tliicknesses shown plotted in Figures 21 to 27. 

Measured Values of Film Thickness 

Meas\ired data is contained in Table III. Pressure readings were 

plotted and the straight line slopes taken to give dp/dz. Frccx Equa- 

tion (71), by neglecting ccropared to p, can readily be computed. 

Similarly can be found from its definition end the measxired values 

of X . Calculations are also included in l^ble 111. 
o 
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Liqxild « 6^1 

Water Orifice: 0.^ in. Jg. 
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